Dynamics of k-core percolation 
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In many network applications nodes are stable provided they have at least k neighbors, and a net- 
work of k-stable nodes is called a k-core. The vulnerability to random attack is characterized by the 
size of culling avalanches which occur after a randomly chosen k-core node is removed. Simulations 
of lattices in two, three and four dimensions, as well as small world networks, indicate that power- 
law avalanches occur in first order k-core systems, while truncated avalanches are characteristic of 
second order cases. 
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The k-core of a network is the set of nodes and edges 
which remain after all nodes which have less than k neigh- 
bors, and the edges attached to them, have been culled or 
removed. Percolation of k-cores has applications ranging 
from magnetism (where it was invented and called boot- 
strap percolation) [il, 01 and rigidity percolation [3, 01, to 
social networks P and protein networks [g^, as well as the 
jamming transition[7,]. The behavior on Bethe lattices 
is particulary interesting due to the mixed nature of the 
transition where a first order jump in the order parameter 
is co-existent with a square root singularity p], This 
behavior has recently been associated with the jamming 
transition in granular media Q and has been analysed 
using a expansion, which confirmed that it persists 
in finite dimensions and by large scale simulations, 
which suggest that in four dimensions the transition is 
not mixed Q. In general, k-core percolation in finite di- 
mensions is quite different than that occuring on random 
graphs, with metastability being a key new feature, so 
that for k > z/2 hypercubic lattices have thresholds ap- 
proaching one in the large lattice limit p^. fill. [T^. [l4| . 
whereas on Bethe lattices the k-core threshold occurs at 
values less than one for all A: < z — Moreover on 

random graphs k-core percolation is first order for all 
A; > 2[1] whereas on triangular lattices and cubic lattices 
the transition is second order and in the connectivity per- 
colation university class for k = sfH, [H, [O, d, fT9l[20|. 

The vulnerability of networks to random node removal 
or attack, has a variety physical applications, including: 
radiation damage; species extinctions in biological net- 
works; and random outages or cascading failures in com- 
munication or transportation networks [2 Ij. There are 
also relations to avalanches in materials problems, includ- 
ing fracture and earthauak esl22 , 23 1, and Barkhausen 
noise in magnetic materials [12" 37 iE]. In fact boot- 
strap percolation is related to the weak disorder limit of 
one of the key models of disordered magnets, the random 
field Ising model (RFIM)!!^]. 

In an earlier paper [l^], we defined elementary culling 
avalanches (EGA) in k-core or bootstrap problems to be 



avalanches induced by random removal of one node on 
a stable k-core, which is a dynamical procedure analo- 
gous to that used in crackling noise simulations in ran- 
dom field systems [25]. We showed that on triangular 
and cubic lattices the EGA distribution provides a sen- 
sitive indicator of the nature of the bootstrap transi- 
tion, even in the slowly convergent second order case 
k = 3 Jl^JM [ll [M [II [13 and the metastable k = 4 

inii,ii3,iii- 
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case 

Here we show that in systems ranging from small world 
networks to regular lattices, including hypercubic lattices 
in four dimensions, the dynamics of k-core percolation 
falls into two classes: (i) Second order systems which 
have truncated avalanche distributions where the largest 
avalanche grows very weakly with system size and; (ii) 
First order systems where avalanches are power law dis- 
tributed and the largest avalanche is of order the system 
size, producing a first order jump in the infinite cluster 
probability at percolation. In the case of Bethe lattices, 
power law avalanches in first order cases can be demon- 
strated explicitly (see below). 

Gulling cascades for a variety of lattices in finite di- 
mensions are presented in Figs. 1 and 2a. These distri- 
butions are found by starting with a lattice with all sites 
present and then repeating two steps: (i) Randomly re- 
move a site and all edges connected to it, (ii) Recursively 
cull all unstable sites until a stable k-core is achieved. 
The number of sites removed during step (ii) constitutes 
an elementary culling avalanche (EGA). The distribution 
of such avalanches, beginning with an undiluted lattice 
and continuing until the lattice is empty yields the cum- 
mulative culling avalanche distribution, C{a). We also 
monitor the largest avalanche that occurs for a given lat- 
tice, Umax, which is the maximum cascade over the whole 
avalanche trace. The finite size scaling behavior of the 
maximum cascade is presented in the insets to Figs. 1 
and 2a. 

In Fig. 1 we present avalanche distributions for k = 
3, 4, 5, on eight-coordinated lattices in three (Fig. la) and 
four dimensions (Fig. lb). These lattices where chosen to 
demonstrate the effect of spatial dimension on the nature 
of the k-core transition. It is clear from this data that 
the fc = 4 transition changes from second order to first 
order as we go from three to four dimensions as noted 
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FIG. 1: (Color online) Efltect of lattice dimension (d), at 
fixed lattice co-ordination (8), on the avalanche distribution 
(EGA), C(a), for the cases k = 3(black circles), k = 4(red 
squares) and k = 5 (blue diamonds), a) Body-centered cubic 
lattices (d = 3) and, b) Simple cubic lattices {d — 4). The 
4-core transition changes from second order for d = 3 to first 
order for d = 4 as is clearly seen in the avalanche distri- 
butions and in the finite size scaling of the largest avalanche, 
a-max (Insets). In a) 1000 samples were used, and in the main 
figure the lattice size L — 128. In b) 1000 samples were used, 
and in the main figure L — 32. 



recently 0, however the k — 3 transition is second order 
in both three and four dimensions. In contrast on z = 8 
Bcthc lattices and random graphs, k-core percolation is 
first order for aU fc > 2. Clearly dimension plays an 
important role on the nature of the k-core transition, 
moreover whenever the transition is first order we see 
the same characteristic behavior: power law avalanches 
and a largest avalanche of order the sample sizep7|. 

We investigated a second method for changing the na- 
ture of k-core transitions, namely by changing the range 
of the interactions. To do this, we created small world 
networks by starting with an undiluted triangular lattice 
and randomly removing a fraction / of bonds from it. 
These bonds are then randomly put back into the lat- 
tice, but they may connect any pair of sites at arbitrary 
distances, with the constraint that all sites remain six- 
coordinated. In this way we create six co-ordinated small 



FIG. 2: (Golor online) k-core percolation on six co-ordinated 
small world networks, with fraction / of small world con- 
nections, a) Avalanche distributions, C(a), for 3-core cases 
with / = 0.15 (blue diamonds), / — 0.20 (red squares) and 
/ — 0.25 (black circles). Inset: The largest avalanche, amax, 
starting at the bottom, for / = 0.15, 0.17, 0.25. b) The infi- 
nite 4-core probability, P™/ versus site probabiltiy p for small 
world fractions (from right to left) / = 0,0.15,0.20,0.25,1. 
Inset: the bootstrap percolation threshold, Pc(/). There is a 
singular behavior in Pc(f) as / — » (see text). In the main 
figures of a) and b) and in the inset to a) from 1000 sam- 
ples for each /, with L — 1028. In the inset to b) from 100 
samples for L = 256(black circles), L = 512 (red squares), 
L — 1024(blue diamonds). 



world networks with fraction / of long range connections. 
When / is zero, we have a regular triangular lattice, while 
when / = 1, we generate a random regular graph with 
co-ordination six everywhere. By varying /, we monitor 
the effect of long-range connections on the nature of the 
k-core transition, which also measures the vulnerability 
of the small world networks to random attack 21]. For 
A: = 3, we find that, as illustrated in Fig. 2a, the tran- 
sition changes from second order to first order when the 
fraction of small world connections reaches fc — 0.23(2). 
For f > fc there are power law avalanches and an ex- 
tensive largest avalanche, while for f < fc, the avalanche 
distribution is truncated, with no large avalanches. 

The 4-core transition on 6-coordinated networks is 
first order for all /, however on triangular lattices it is 
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metastable so that for large lattices the threshold ap- 
proaches one very slowly, while on random graphs there 
is a finite threshold Pc < I- Nevertheless we find power 
law avalanches and an extensive largest avalanche in all of 
these cases. However, the behavior of the k-core thresh- 
old as a function of the fraction of small world connec- 
tions is interesting and is presented in Fig. 2b. A no- 
table feature of this data is that small world connections 
quickly stabilize the metastable 4-core of triangular lat- 
tices and, as is emphasized in the inset to Fig. 2b, there is 
a singular behavior in pdf) as / — > 0. Presumably, this 
non-perturbative behavior is related to the stabilization 
of large holes in the triangular lattice by the addition of 
long-range loops characteristic of small world networks, 
and it may be possible to characterize this behavior using 
rigorous mathematical methods [l^, [ll|, [l^, [l^ . 

We also observed power law avalanches and an exten- 
sive largest avalanche in the limit / = 1 which is the 
random graph limit. In this limit power law avalanches 
can be demonstrated explicitly using Bethe lattices, as 
we now show. Consider site diluted, z-coodinated Bethe 
lattices where a fraction p of the nodes are randomly re- 
moved, and let T be the branch probability that a node 
of the Bethe lattice is part of a k-core. The branch proba- 
bility, T, is found by solving the self-consistent equation, 

T = pY, K-'T-- wtth A\ = (1 - Tf-' (1) 

n=k ^ ^ 

which is well known from the original bootstrap paper Q 
and from applications to rigidity percolation The 
branch probabilities are used to find the probability of 
a site being on the k-core, 

z 

X^pY^AIT-. (2) 

n—k 

We build on the results (1) and (2) to find the culling 
avalanche distributions for a branch, B(a), and for the 
Bethe lattice, P{a). P{a) is the differential probabil- 
ity and is related to the cummulative probability, C(a) 
through C{a) = / P{a,p}dp^_^ 

As noted previously |28|, |29|, sites which have exactly 
k neighbors play a special role in culling avalanches on 
Bethe lattices. Following their terminology, we define 
the corona as a connected cluster of sites each of which 
is part of the k-core and has exactly k neighbors. With 
this definition, removal of a site of the k-core leads to 
a culling avalanche in the corona surrounding that site. 
The probability of finding a corona of size, a, surrounding 
a randomly chosen starting site on the Bethe lattice is 
given by, 

z n 

P{a)=pJ2All[Y.B{ai)S{S,,). (3) 

n—k l — l ai 

where Sn = Ym=i — a + 1. This expression has the 
same form as Eq. (2), but with the term T" in that 
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FIG. 3: A log-log plot of culling avalanches, P{a) for 3-core 
percolation on a z—4 Bethe lattice, found from P*{y) (see 
text after Eq. (8) using an exact Mathematica algorithm. 
The dark lines are for 5 values of the site concentration, from 
the top starting with p = 8/9 = pc and then moving above pc 
in increments of Sp = 10~*. The continuous line of slope 3/2, 
vertically shifted from the data, is included to demonstrate 
that P(a) ~ 1/a^^^ at the 3-core threshold, pc- 

equation replaced by the product in Eq. (3). This term 
is a product over the probabilities of finding corona on 
each of the branches connected to the starting site, and 
the delta function ensures that the ECA's on these corona 
sum to the total avalanche size a. From this equation 
it is evident that the normalizations ^ Pio.) = X and 
^i?(a) = T hold. A similar reasoning is applied to the 
branch probabilities, yielding the extension of Eq. (1) to, 

* = ^n-'T-Sia) + Al-_\ n E B{a,)6{S,-,). 

^ n=k 1=1 ai 

(4) 

There are two terms in this expression, as only corona 
sites on a branch initiate a culling avalanche on the 
branch. All other k-core sites terminate the EGA and 
contribute to the prefactor of S{a) in this equation. This 
delta function is not present in Eq. (3) as the starting site 
there is removed by hand to initiate an EGA so it can 
have any co-ordination. Defining the generating func- 
tions, B{y) = E^Wy", and P{y) = and 
applying them to Eq. (3) yields, 

z 

P(y)^PvY^AlB{yT. (5) 

n— 

Similarly, Eq. (4) becomes, 

z-l 

B{y)^pY^Al-^T^^pyAl7_\B{yf-\ (6) 

n— /c 

The probabilities B[a) or P{a) are extracted from the 
generating functions B(y) or P{y), by differentiation or 
by contour integration, whichever is more convenient. 

The cases fc = 1, 2 are the standard connectivity perco- 
lation and percolative backbone problems respectively, so 
the simplest non-trivial k-core problem on random graphs 
is fc = 3, z = 4. The formalism outlined above can be 
carried through explicitly in this case and has a behavior 
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characteristic of all cases > 3 [2 81 129[. For this case, 
B{y) obeys, 



B{y) = pT^ + 3py{l - T)B{yY 



(7) 



where T = 3/4 + 9((5p)^ /8a/2, which follows from Eq. (1) 
and Sp = p — Pc- This illustrates the characteristic fea- 
ture of a jump discontinuity and a continuous square root 
singularity, as noted before[ll,Q. Solving Eq. (7) and se- 
lecting the physical root, we find an explicit expression 
for the avalanche generating function, 



1 - [1 - Up^T^il ~ T)y]i 
6p(l - T)y 



(8) 



and finally, from Eq. (4) , we find the generating function 
for Bethe lattices, P*(y) ^ pyl^il - T)B*{yf + B* (y)^]. 
We wrote a Mathematica program which finds the exact 
culling avalanche distribution from P*(y), and the re- 
sults are presented in Fig. 3, which demonstrates the 
3/2 power law behavior at threshold characteristic of 
mean field avalanches. Asymptotic analysis near crit- 
icality is carried out by noting that the normalization 
condition B*{y = 1) = T, enables an expansion of B* 
near y = 1 and p — Pc- The leading term in the expan- 
sion of B* is then, B* 1 - [1 - y{l + 96p/2)]^/'^, from 
which we find that near the threshold, Pc, the tail of the 
avalanche distribution obeys P(a) ~ [l—D{p—pc)]"'/a'^^^, 



where D = 9/2 and Pc = 8/9 for 3-core percolation 
on 4-coordinated Bethe lattices. In experiments as well 
as computer simulations, it is more convenient to mea- 
sure the integrated probability distribution of avalanches 
C(a) and in the simulation results presented above, we 
integrate over the complete avalanche trace. The proba- 
bility of large avalanches is significant in a small region 
of width near threshold. Integration of P{a) in this 

regime yields, C(a) ~ a^^, for a — > oo which is ob- 
served in our simulations on random graphs, and is typ- 
ical of avalanches occuring in infinite range models and 
Bethe lattices in a wide variety of non-equilibrium sys- 
tems f2i,[2l[2i,[3S[3l. 



In all cases we studied, power law avalanches are 
associated with first order k-core transitions, both in 
metastable and regular first order cases. This is not 
without precident as in random field Ising hysteresis a 
first order jump occurs in the hysteresis loop despite the 
fact that power law avalanches occur [l^, How- 
ever in brittle fracture problems in finite dimensions, 
the avalanches are truncated when the loading sharing 
is short ranged [33. [33 |. Presumably this is due to the 
greater amplification of stress on survivors in fracture 
problems with short range load sharing. Nevertheless, 
k-core dynamics capture the essential physics of instabil- 
ities in a broad class of complex network problems. 
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